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bo"' Abstract 

O. 

We discuss the open charm production in peripheral reactions pp — » Y C Y C and pp — » M C M C , 
where Y c and M c stand for A+,X+ and D,D*, respectively, at -^/i < 15 GeV, which corresponds 
to the energy range of FAIR. Our consideration is based on the topological decomposition of 
the planar quark and diquark diagrams which allows to estimate consistently meson and baryon 
exchange trajectories and energy scale parameters as well. The spin dependance is determined 
by the effective interaction of lowest exchanged resonance. Unknown parameters are fixed by an 
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independent analysis of open strangeness production in pp — > YY and pp —> KK reactions and 
of SU(4) symmetry. We present the corresponding cross sections and longitudinal double-spin 



asymmetries for exclusive binary reactions with open charm mesons and baryons in the final state. 
The polarization observables have a non-trivial t and s dependence which is sensitive to details of 



the open charm production mechanism. 
PACS numbers: 12.38Bx, 13.85M, 13.88. +e 



I. INTRODUCTION 

Open charm production will be one of the major topics of the hadron and heavy-ion 
programme at FAIR [1]. On the one hand, charm spectroscopy will be addressed by the 
PANDA collaboration [2], while the CBM collaboration [3] will exploit charmed particles as 
probes of the nuclear medium at maximum compression. For both large-scale experiments 
at FAIR one needs to know the properties of charmed baryons and mesons as well as their 
production processes in elementary pp and pp reactions. For this the opportunities at FAIR 
are promising, as for instance, the PAX collaboration [4]) envisages the use a polarized 
anti-proton beam. This offers the chance to study in depth the mechanism of open charm 
production at the moderate energies from threshold to ^/s < 15GeV. In this energy range 
the phenomenology of charm production is not well established. In present paper we select 
one important problem of this wide field, namely, the analysis of exclusive binary reactions 
pp — ► Y C Y C , pp — ► DD, pp — ► DD* etc., in peripheral collisions in the mentioned energy 
range. 

Since the initial energy is not asymptotically high, the widely used models for the heavy 
quark production based on perturbative QCD (see for example Ref. [5-7]) are not applicable, 
and an essential improvement by including high order corrections is needed [8, 9]. Another 
severe problem is related to the dynamics of charm productions. In the popular QCD 
models, the c quark is produced through gluon fragmentation. For c production in peripheral 
collisions such a gluon must have a large momentum (large x ~ 1), i.e., much larger than its 
average value in a nucleon x < 0.2, and therefore, this mechanism is strongly suppressed. 

It is expected that the mechanisms of peripheral charm production in pp — ► DD and 
pp — ► Y C Y C reactions are similar to the strangeness production in pp — > KK and pp — ► YY 
reactions, respectively, which were described successfully in terms of Regge pole models [10- 
13] with certain baryon and vector meson exchange trajectories. However, a direct extrapo- 
lation of such models to the charm production is faced by a problem. First, the linear Regge 
trajectories leads to large negative values of the corresponding intercepts «a c (0) ~ —4.5 and 
Od*(0) ~ —2. This would result in a suppression of the charm production in peripheral 
reactions which contradicts to the corresponding data on inclusive charm production. This 
means, in turn, that the trajectories connected to masses and spins of the charmed hadrons 
must be essentially non-linear (cf. Refs. [14, 15]. Another problem is the estimate of the 



energy scale parameter in the Regge pole propagator, which also affects much the cross 
sections. 

Therefore, it seems to be interesting and important to use a model based essentially on 
a non-perturbative QCD background being reliable for describing the peripheral reactions. 
Such an approach was developed in Refs. [18-20] and applied for the evaluation of cross 
sections of the exclusive A c production in np and pp collisions. The binary ir~p — > D~A C 
exclusive process plays an important role in this consideration [19]. The model for this 
reaction is based on quark-gluon string dynamics, assuming the annihilation of a qq pair 
in the interaction, the formation of a qq color tube with subsequent decay to the observed 
hadrons (see also Refs. [21, 22]). Schematically, the process n~p — ► D~A C is described by 
the planar diagram shown in Fig. 1 (a). A more general case is exhibited in Fig. 1 (b). 
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FIG. 1: Diagrammatic representation of planar diagram for reaction n~p — ► D~A C (a) and for 
more general case (b). 

The assumption of the formation and decay of color-gluon strings allows to construct the 
space-time evolution of the process and to obtain the factorization condition, where the 
imaginary part of the amplitude of the process ab — > cd is expressed via a product of the 
probabilities w a b and w c d of the elastic scattering of ab — > ab and cd — > cd, respectively. This 
gives a consistent prescription for evaluating the parameters of the amplitude (trajectories 
and energy scale parameters) of the non-diagonal transition ab — »■ cd. 

The aim of the present paper is to extend the results of Refs. [18-20] for exclusive charm 
production in the pp collisions. In our consideration we analyze simultaneously the open 
charm and open strangeness production. We are going to consider the reactions pp — » AA 
and pp — > A C A C , pp — > KK and pp — ► DD etc. The strangeness production has its own 
interest, but on the other hand, some available (although relatively old) experimental data 
allow to fix the unknown parameters of the model and get absolute values of the cross sections 
of the open charm production. We also analyze the double longitudinal asymmetry which 
will be accessible in mentioned the FAIR PAX experiment due to the polarized anti-proton 



beam. The spin dependence of the amplitudes is generated by the symmetry of the PNY 
and VNY interactions (P — K,D and V = K*,D*...) which was widely used in various 
studies (see e.g. Refs. [23, 24]). 

For completeness, we mention that some aspects of the inclusive charm production in 
terms of quark-gluon string model were discussed in Refs. [25, 26]; polarization effects in 
open charm photo-production were considered in Ref . [27] ; propagation of charmed hadrons 
in the nuclear medium were analyzed in Refs. [28-30] ; open charm production in relativistic 
nucleus- nucleus collisions at wide energy region was analyzed in Ref. [31]. 

Our paper is organized as follows. In Sec. II we analyze the strangeness production in 
the reactions pp — ► AA, pp — > AX°, and pp — »■ X°X°, and the open charm production in 
pp — ► A+A+, pp — ► A+X+, and pp — ► X+X+ processes, where the dominant contribution 
comes from the K* and D* exchange trajectories, respectively. First, we note equations for 
the invariant amplitudes and then discuss our results for the differential cross sections and 
the longitudinal asymmetry. In Sec. Ill we provide a similar analysis for the pp — > KK 
and pp — ► DD reactions, assuming the dominance of strange and charmed baryon exchange 
trajectories. The reactions pp — > KK* and pp — > DD* are discussed in Sec. IV. The 
summary is given in Sec. V. 

II. REACTIONS pp -► YY AND pp -► Y C Y' C 

In this section, we discuss strange and charmed baryon-antibaryon production in periph- 
eral pp collisions. For the sake of simplicity, we consider the exclusive production of AA and 
A C A C pairs. The generalization for reactions with AE, SA, XX final states may be done in 
a straightforward manner. 

The corresponding planar diagrams for AA and A C A C are depicted in Fig. 2 (a) and (b). 

^ d_ ^d t - d_ ^d 7 

P u i la A P a I l - n A c 

I) C . c 



u A D U » -y- U A 



d- -d v d- ^d 



FIG. 2: Planar diagram for the reactions p — ► AA (a) and p — > A C A C (b). 



A. Reaction pp — ► AA 

Following Ref. [19] we assume that the amplitude of the reaction pp — > AA has the form 
of a Regge pole amplitude, dominated by the K* exchange trajectory, 

g 2 » / s \ asq ^~ 1 

s \ s pp:AA J 

where mi, m/, rii and n/ are the spin projections of p, A, p, and A, respectively, q stands 
for u and d quarks, asq(t) is the K* + trajectory, g 2 K * Nlv is the coupling constant of the 
K*NA interaction, and sq = 1 GeV is an universal scale parameter. The spin dependence is 
accumulated in the amplitude M. which is a smooth function of the Mandelstam variables 
s and t. In the limit of s — ► oo one has M. oc s. The explicit form of M. will be defined 
later on. The overall residual function C(t) will be found from a comparison with available 
experimental data. 

In our consideration we use the nonlinear representation for the meson trajectories de- 
veloped in Ref. [15] 

«(t) = a(0)+ 7 (v / T-v / T^7), (2) 

where 7 = 3.65 GeV -1 is the universal parameter (a universal slope in the asymptotic 
region), and T ^> 1 GeV 2 is the scale parameter, being special for each trajectory. In the 
diffractive region with —t <^C T, the linear approximation 

a(t) = a(0) + a't, (3) 

is valid with a' = , y/2yT. 

The intercept a Sq (0) and the slope a' s of the trajectory for the non-diagonal transition 
are related to the corresponding parameters for diagonal transitions as following [15, 18] 

2a- sq (0) = a qq (0) + a s - s (0) , (4) 

2/a'- sq = l/a' qq + l/c4 , (5) 

where (Xq q (t) and a 5s (t) are the p and <fi meson trajectories, respectively. 
In our numerical calculations we employ 



a p (0) = 0.55, t/T~ p = 2.46 GeV, a' ~ 0.742 GeV 



-2 



o^(0) = 0.414, vTk* = 2 - 58 GeV , a K* - 0-71 GeV" 2 , 

cty(O) = 0.28, y% ~ 2.70 GeV, a', ~ 0.676 GeV" 2 , (6) 



where the p and K* trajectories are taken as input according to Ref. [15]. 

The energy scale parameter s pp . AA in Eq. (1) is related to the corresponding scale param- 
eters for the diagonal transitions pp — > pp, (s pp ) and AA — ► AA, (saa) as 

(^ : aa) 2( ^* (0K1) = (s Pp r {0) - 1 x (saa)^- 1 • (7) 

The scale parameter for the diagonal transition s a b is determined by the sum of the transver- 
sal masses of the constituent quarks [19] 

(n a \ / n b \ 

with M q± ~ 0.5 GeV, M s± ~ 0.6 GeV, and M c± ~ 1.6 GeV. This leads to the following 
values for the scale factors: s pp ~ 2.25 GeV 2 , saa — 2.56 GeV 2 , and s Pp .xa — 2.43 GeV 2 . 

We assume that the spin dependence of the amplitude in Eq. (1) is determined by the 
symmetry of the NAK* interaction given by the effective Lagrangian in the conventional 
form 

C K . NY = -Y ( 7 „ - j^f^^) NVIC* + h.c. , (9) 

where N, Y and K* denote the nucleon, hyperon and the K* meson fields, respectively, Y 
stands for A, S, etc., and k is the tensor coupling strength. 
Using this form, one can get the amplitude Ai in Eq. (1) as 

■^m f n f ;minf\ S ^) = ^ {S,t) L mfm . 1 n ^ n . {— Q^ v -\ ~2~) > {*-") 

where q is momentum transfers in the pAK* vertex: q = p p — p\, with p p and p\ as four- 
momenta of the incoming proton and outgoing A. The functions T^ p '^ read 



r ( P ) - fft, \ (Pp+Pa) 

^ = U A \{l + KK'NAhn-KK'NA M +M 
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Tf = v P ((1 + Kk-naH + "x-^ m^+m'J ) VK • (11) 

The normalization factor Af(s, t) eliminates additional s and t dependence provided by the 
Dirac structure in Eq. (10) which is beyond the Regge parametrization: 

A/"M) = |^|, F oc (s) = 2s, 

F 2 (s,t) = Tr (r^T^"'*) Tr (r^T®" 4 ) (g^ - ^)(<V*' - Qj ^f) . (12) 



For the NY K* coupling constants we use the average values of the Nijmegen poten- 
tial [32]: g K * NY = -5.18, k k *ny = 2.79 for Y = A and (-3.29, -0.91) for Y = E°. 
The cross section is related to the invariant amplitude of Eq. (1) as 

da 1 



dt 16tt(s - AM 2 N ) 2 



T fi \ 2 , (13) 



where summing and averaging over the spin projection in initial and the final state is pro- 
vided. We will also discuss the longitudinal double spin asymmetry, defined as 

do* - da^ 
A da^ + da*' [ } 

where the symbols *=> and =4 correspond to the anti-parallel and parallel spin projections of 

incoming p and p with respect to the quantization axis chosen along the proton momentum 

in the center-of-mass system (c.m.s.). 

The generalization for reactions pp — > AE, pp — > EA, and pp — ► EE is accomplished by 

the substitution of Ma — >■ Ms, Qk*nk -^ 9k* ivs and kk*na — > k*k*nt, in Eqs. (1) and (11). 

B. Reaction pp — ► A C A C 

In this case, the amplitude is defined by Eq. (1) with the obvious substitution A — ► A+ = 
A c , E° — >■ E+ = E c , K* — >■ D*, a,^ — > oljj^, etc. As a first approximation, we assume the 
validity of SU(4) symmetry and, therefore, the coupling constants of the D*NY C interaction 
are chosen to be the same as for the case of K*NY interaction. The corresponding trajectory 
and the energy scale parameters read 



a D *{0) 


= -1.02, 


y/T D . = 3.91 GeV, 


ttj/V>(0) 


= -2.60, 


y/T J/lP ~ 5.36 GeV 


S pp:A c A c 


~ 5.98 GeV 2 . 
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a' n , ~ 0.467 GeV" 



i D * 



a'jfr ~ 0.34 GeV" 2 



(15) 



1. Differential cross sections 

Consider first the strange hyperon production pp — ■» VV which we use to fix the residual 
factor C(t) in Eq. (1). In Fig. 3 we show the differential cross section of the reaction 



pp — > A A and pp — ► AS as a function of the momentum transfer £ = (p p — Py) 2 at the 
initial momentum pl = 6 GcV/c together with the available experimental data [38]. The 
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FIG. 3: Differential cross section of 
the reactions pp — ► A (solid curve) and 
pp — ► AS (dashed curve) as a func- 
tion of the momentum transfer t at 
PL = 6 GeV. The experimental data are 
taken from Ref. [38]. 
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provides a reasonable agreement of the calculation and the data. 

In Fig. 4 (left panel) we exhibit our prediction for the differential cross sections of the 
reactions pp — > AA, pp — > AS (E°A) and pp — > E°E° as a function of t max — t at initial 
momentum p^ = 10 GeV/c. Here, t max is the maximum momentum transfers which corre- 
sponds to the A production at zero angle relative to the momentum of the incoming proton 
in the c.m.s. 

The exponential decrease of the cross section is defined by the Regge propagator 
(s/si) 2aK *^ and the residual C(t). The dependence on the excess energy As 1 / 2 = \^s — ^/s^, 
where ^/sq = My' + M Y , is shown in Fig. 4 (right panel). The calculation is done at fixed 
tmax — t = 0.2 GeV 2 . At large energies, the cross section behaves as s 2 *"^* -1 ) ~ s -1 - 172 . The 
ratio of the cross sections with AA, AE° and E°E° final states at large energy reads 



1 : r : r A 



(17) 



where r = (g K * NA / 'g K *NY.) 2 ^ 0.4. 

The predicted differential cross sections of the charm hyperon production are exhibited 
in Fig. 5. Here, we use the notation A c = A+ and E c = E+. The left panel shows the 
dependence on t max — t at fixed pi = 15 GeV/c. The right panel exhibits dependence on 
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FIG. 4: Left panel: The differential cross sections of the reactions pp — ► AA (solid curve), pp — > AS 
(dashed curve), and pp — ► S°S° (dot-dashed curve) as a function of t max — t at pl = 10 GeV/c. 
Right panel: The differential cross section as a function of the excess energy As 1 ' 2 at t msx — t = 
0.2 GeV 2 . 
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FIG. 5: Left panel: The differential cross sections of the reactions pp — ► A C A C (solid curve), 
pp — ► A C S C (dashed curve), and pp — ► S c S c (dot-dashed curve) as a function of t max — t at pi, = 
15 GeV/c. Right panel: The differential cross section as a function of the excess energy As 1 ' 2 at 
W-t = 0.2GeV 2 . 



the energy excess As 1 / 2 at fixed t max — t = 0.2 GeV 2 . The threshold initial momenta (for 
a fixed target) for reactions with A C A C , A C X C and S C S C final states are 10.15, 11.83, and 
10.85 (GeV/c), respectively. The energy excess at pi = 15 GeV/c for these final states are 
0.90, 0.571 and 0.570 GeV, respectively. This energy is not asymptotically high and some 



particular behavior of the cross sections in the pre-asymptotical region is expected. Thus, 
in Fig. 5 (right panel) one can see a bump-like behavior at low As, which reflects the energy 
dependence of t max in this region. 

2. Longitudinal asymmetries 

For a better understanding of the results of our numerical calculation, it seems to be useful 
to perform a qualitative analysis of the longitudinal asymmetry at forward production angle 
(or t = t max ), where the orbital interaction is absent. In this case, the amplitude of the 
pp — ► YY reaction may be written as 



( A(s) 5 mimf 5 ninf + —=B{s) (1 - 477ljTO/) b~~„ Hmf O~-n ins J 



T r n s n r ,m t ,n t ~ R(s) ( A(s) 5 mimf S ni n f + —=B{s) (1 - 477ljTO/) 8- mimf 8 

where R(s) is a spin-independent function, and mj,m/, 77j and 77/ stand for the spin projec- 
tions of p, Y, p, and Y, respectively, The longitudinal asymmetry is expressed through the 
spin-conserving (A(s)) and spin-flip (B(s)) amplitudes as 

A^a) + B 2 (s) ' ( ' 

The spin- conserving amplitude is determined by the two functions a (s) and a K (s) 

A(s) = (oo(s) + a K (s)) 2 , 
1 PpPy 



a (s 



(E + M N )(E + M Y ) ' 



f \ f \ 2kE (i p P py ^ (on\ 

a K (S) = KCLnls) 1 — ^- , 20 

kK ' 0K } M N + M Y \ (E + M N )(E + M Y )J K ' 

where k is again the tensor coupling strength; p p , py denote the three momenta of the 
proton and the outgoing hyperon, respectively, and E = \fs/2 is the proton energy in c.m.s. 
In the case when M p ~ My or/and at high energies, when yfs ^> My, a K — >■ and the 
spin-conserving amplitude becomes independent of k. 

Contrarily, the spin-flip amplitude is proportional to the square of the magnetic strength 
(1 + k) 2 : 

B w=-^(< i+ «><Ffk-^>y '■ ' 2i > 

At high energies with E ^> My and p p ~ py, B(s) — > and, therefore, the asymmetry in 
Eq. (19) vanishes. However, at finite energies and large values of (1 + k), the amplitudes 
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FIG. 6: The quantities A 2 (spin-conserving amplitudes, solid curves) and B 2 (spin-flip amplitudes, 
dashed curves) as a function of the excess energy As 1 ' 2 with t = t max at different values of the 
tensor coupling k. The left and right panels correspond to the reactions pp — ► AA and pp — ► A C A C , 
respectively. 

A(s) and B(s) are comparable, and the longitudinal asymmetry may be finite and large. 
The energy dependence of A 2 (s) and B 2 (s) for the reactions pp — >■ AA and pp — * A C A C is 
shown in Fig. 6. In the AA final state, the function B 2 (s) is rather small due to the small 
difference between Mn and M\. The dependence of A 2 (s) on the tensor coupling k is rather 
weak. This leads to the small value of the longitudinal asymmetry for the pp — > AA reaction. 

In case of the pp — > A C A C reaction, the situation is quite different. The large difference 
between M^ and M\ c leads to a large value of B 2 (s), shown in Fig. 6 (right panel), and 
results in a large value of the longitudinal asymmetry. 

For the pp — ► AS (SS) reactions the spin-flip amplitude |-B(s)| is small because of the 
small magnetic strength, 1 + k ~ 0.09, and the asymmetry is almost zero. 

Our predictions for the pp — > YY reactions are shown in Fig. 7. The left panel exhibits 
the t dependence at the initial momentum pi = 10 GeV/c. The right panel shows the 
dependence on the energy excess at t max — t = 0.2 GeV 2 . One can see that the result of 
the numerical calculations agrees with our qualitative consideration. Thus, for pp — > AA, 
the asymmetry does not exceed 0.2 at forward angles and decreases with energy. In the 
pp — > AS and SS reactions it is almost zero. 

The longitudinal asymmetry for the pp — ► Y C Y C reactions is presented in Fig. 8. In the 
case of the pp — > A C A C reaction, the asymmetry is large at low energy excess and decreases 
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FIG. 7: Left panel: The longitudinal asymmetry for the reactions pp — ► AA, AS and S°S° as a 
function of momentum transfer £ at j>l = 10 GeV. Right panel: The asymmetry as a function of 
the excess energy at t max — t = 0.2 GeV 2 . 
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FIG. 8: Left panel: The longitudinal asymmetry for the reactions pp — ► AA, AS and A°S° as a 
function of momentum transfer t at pl = 10 GeV. Right panel: The asymmetry as a function of 
the excess energy at t max — t = 0.2 GeV 2 . 

rapidly with energy. In the reactions pp — > A C S C and S C S C the asymmetry is negligibly 
small. 
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III. REACTION pp -► MM 

In this section, we discuss the production of MM (with MM being AA or DD) in 
pp collisions. We assume that at small momentum transfer —t, where t = (p p — Pk) 2 or 
^ = iPp ~ Pb) 2 i the dominant contribution comes from the baryon exchange channels. 

As an example, in Fig. 9 (a) and (b) we show the planar diagrams for pp — > A~A + and 
pp — ► D°D°, with A and A+ exchange, respectively. The cases of E (S+) exchange, or E + 
(E+ + ) exchange for AA (D + D~) final state, are similar. Here and further on we employ 
the quark-diquark identity, used in many phenomenological approaches to QCD [14, 39, 40]. 
This means that the exchanged baryon is considered as a quark-diquark string object, shown 
schematically in Fig. 9 (c). 

K(D) 



K(D) 



FIG. 9: Planar diagrams for the reaction pp — ► K K + (a) and pp — ► D°D° (b). The exchanged 
baryon as a quark-diquark object (c). The symbol d stands for a qq diquark. 



A. Reaction pp — ► KK 

Let us consider first the reaction pp — > KK. For definiteness, we consider pp — > A~A + 
with A exchange. The cases of E and E + for K°K° final states may be executed in an analog 
way. 

The assumption of the quark-diquark identity allows to generalize our model developed 
in the previous section. Namely, we assume that the amplitude of the reaction pp — > K~K + 
has the form of a Regge pole amplitude dominated by the A + E exchange trajectories. Thus, 
for the A exchange it reads 

TZ-?- K+ = C\t)MZ^ K \s,t) 9 ^T^-a d M) — , (22) 

s £ \ s pp:KK J 

where rrii and nj are the spin projections of p and p, respectively, d stands for a ud diquark, 
cxds(t) is the A trajectory, g 2 KNlv is the coupling constant of the AAA interaction and s = 

13 



1 GeV is an universal scale parameter. The spin dependence is accumulated in the amplitude 
A4 which in the limit of s — ► oo results in M. ?a ^fs. The explicit form of M. will be 
defined later. The overall residual function C'(t) will be found again from a comparison 
with available experimental data. 

The parameters of the trajectory for the non-diagonal transition a ds are related to the 
corresponding parameters for the "diagonal" transitions a Ss and a dd similarly to Eqs. (4) 
and (5) 

2a ds (0) = a dd (0) + a s - s (0) , (23) 

2/a' ds = l/<4 + l/a'- ss . (24) 

Using the A trajectory as input [10] 

a ds = a A = -0.65 + 0.94t , (25) 

and a Ss from Eq. (6), one can evaluate the diagonal a dd trajectory at small \t\ as 

a dd (t) = -1.58 + a' dd t (26) 

with a' M = 1.542 GeV" 2 . 

The equation for the energy scale parameter Spp.jcK is slightly different from Eq. (7). Now 
it reads 

/ _ \2(a ds (0)-|) _ / \a dd {0) ( \a Ss (0)-l / 9? x 

\ s pp:KK) — K^pp) x \ b KK) ■ \ ZI ) 

Using Sk K = 1.21 GeV 2 and s Pp = 2.25 GeV 2 , one gets s Pp .ji K = 1.853 GeV 2 . 

The spin dependence of the amplitude in Eq. (22) is determined by the form of the KNK 
interaction given by the effective Lagrangian 

Cnyk = -iN 75 YK + h.c. , (28) 

where N, Y and K denote the nucleon, hyperon and the K meson fields, respectively, Y 
stands for A, E, and so on. This form leads to the following expression of the amplitude M. 
in Eq. (22) 

M^ RK (S, t) = Af( S , t) [V ni {i>y - My) U mi ] , 

A/XM) = f^|, Fl(s) = sM Y /2, 

F2 (s,t) = ^ ((s- 2M 2 N )(M Y -t) + AM N M Y {M 2 N + M 2 K + t 
-(M 2 N -M 2 K + t) 2 -M 2 N (M Y + t)) , (29) 
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where py = p P — Pk- 

For the KNY coupling constants we use the average values of the Nijmegen potential [32]: 
g KNY = -15.755, for Y = A and -4.785 for Y = E°. 

In case of the pp — ► K~K + reaction, the total amplitude is the coherent sum of the A and 
E exchange trajectories, while the reaction pp — > K°K° is dominated by the E + trajectory. 
Following [10] we use 

a s ~ a s + ~ -0.79 + 0.87*. (30) 

For simplicity, for the E exchange channels we use the same the energy scale parameter as 
for the A exchange, taking into account similarity of the corresponding trajectories and the 
fact that the contribution of E exchange amplitude is much smaller than the dominant A 
exchange one. 

B. Reaction pp —>■ DD 

In this case, the amplitude is defined by Eq. (22) with the substitutions A — ► A+ = A c , 
E° -► E+, E+ -h. E++ K+ -h. D°, K~ -+ D°, K° -+ D~ and K° -+ D+ and so on. As above, 
we assume the validity of SU(4) symmetry which means that the coupling constants of the 
DNY C interaction are chosen to be the same as for the case of KNY interaction. The A c 
trajectory is calculated using 

2a dc (0) = a dd (0) + a a (0) , (31) 

2/a' dc = l/c4 + l/o4 , (32) 

where acdt) = eej/^(t) and (Xdd(t) are defined by Eqs. (15) and (26), respectively. Thus, for 
a<ic(t) and the energy scale parameter s Pp , D f) we have 

o4(0) ~ -2.09, a' dc ~ 0.557 GeV~ 2 , (33) 

s pp:DD ~ 3.59 GeV 2 . (34) 

For simplicity, we assume a A + = a dc — a s + — a s ++. 
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C. Results 



1. Differential cross sections 



The differential cross section of the pp — > K~ K + reaction as a function of the momen- 
tum transfer t = (p p — Pk+) 2 at initial momentum pl = 5 GeV/c together with available 
experimental data [41] is presented in Fig. 10. The separate contributions from A and S° 
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FIG. 10: Differential cross section of 
the pp — ► K~K + reaction as a function 
of momentum transfer t at pl = 5 GeV. 
The contributions from A and £ ex- 
changes are shown by dashed and dot 
dashed curves, respectively. The exper- 
imental data are taken from Ref. [41]. 



exchange are shown by dashed and dot dashed curves, respectively. The solid curve is the 
coherent sum of these contributions. One can see a dominance of the A exchange trajectory 
in pp — > K~K + . This reaction is used to fix the residual factor C'(t) in Eq. (22). We find 

0.52 



C'{t) 



(35) 



(l-t/1.15) 2 ' 

i.e., it coincides within ~ 30% with the residual in pp — ► YY reactions (cf. Eq. (16)), which 
is in favor of the consistency of the model. 

In Fig. 11 (left panel), we show our prediction for the differential cross sections of the 
reactions pp — » K~K + , and pp — > K°K° as a function of t max — t at the initial momentum 
Pl = 10 GeV/c. The dependence of the differential cross sections on energy (a/s) at fixed 
t max -i = 0.2 GeV 2 , is exhibited in Fig. 11 (right panel). At large energies the cross sections 



behave as ~ s 



-2.3 



and 



~ s 



-3.58 



for the K K + and K K final states, respectively. 



Our prediction for the differential cross sections of DD pair production is presented in 
Fig. 12. The left panel illustrates the dependence on t max — t at fixed pi = 15 GeV/c. 
The right panel exhibits the dependence on the excess energy As 1 / 2 at fixed t max — t = 
0.2 GeV 2 . The ratio of the cross sections with D~D + and D°D° final states is close to 
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FIG. 11: Left panel: The differential cross sections of the reactions pp — ► K~K + (solid curve) 
and pp — ► K°K> (dashed curve) as a function of t max — t at p^ = 10 GeV/c. Right panel: The 
differential cross sections as a function of the energy y/s at t max — t = 0.2 GeV 2 . 
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FIG. 12: Left panel: The differential cross sections of the reactions pp — ► D°D° (solid curve) and 
pp — ► D~D + (dashed curve), as a function of tmax-t a t PL = 15 GeV/c. Right panel: Differential 
cross section as a function of the excess energy As 1 ' 2 at t max — t = 0.2 GeV 2 . 

(\/2 gKNT,/ 9kna) 4 — 0.034. The cross sections decrease rapidly with energy as s~ 6 - 18 , there- 
fore, the region with small excess energy is more suitable for studying these reactions. 
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2. Longitudinal asymmetry 

In reactions pp — >■ KK {DD) at forward production angle (or t = t max ), the spin in the 
final state is equal to zero. This means that the production amplitude may be expressed as 



J-m.i,rn ~ -t3\ s ) Omi-rii ■> 



(36) 



and therefore, the asymmetry in Eq. (14) A = 1. At finite angles, the spin-orbital interac- 
tions becomes important, which leads to an increase of the contribution of da^ to the total 
cross section and to a decrease of the asymmetry. 
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FIG. 13: The longitudinal asymmetry for reactions the pp — ► K~K + (solid curves) and pp — ► 
K°K°, (dashed curves). Left panel: The asymmetry as a function of momentum transfer £ max — t at 
PL = 10 GeV. Right panel: The asymmetry as a function of the energy y/s at t max — t = 0.2 GeV 2 . 



In Fig. 13 we show our prediction for pp — > KK. The left panel exhibits the t dependence 
at pl = 10 GeV/c, whereas the right panel exhibits the \fs dependence at t max — t = 
0.2 GeV 2 . One can see a decrease of A with —t and an almost constant value at large \/s 
and fixed t max — t. Some difference in A for K°K° and K + K~ final states is explained by the 
difference of masses of K° and K ± mesons, which leads to the different relative momenta 
and to some difference in spin-orbital interactions. 

The longitudinal asymmetries for the pp — > DD reactions are presented in Fig. 14. The 
left panel demonstrates the t dependence at the initial momentum pl = 15 GeV/c. The right 
panel shows the dependence on the excess energy at t max — t — 0.2 GeV 2 . One can see that 
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FIG. 14: The longitudinal asymmetry for pp — ► D°D° (solid curves) and pp — ► D~D + , (dashed 
curves). (Left panel: The asymmetry as a function of momentum transfer t max — t at pi = 15 GeV. 
Right panel: The asymmetry as a function of the energy excess As 1 ' 2 at t max — t = 0.2 GeV 2 

the results of the numerical calculation agrees with our above qualitative consideration. The 
difference in the asymmetries for KK and DD final states is mainly due to the difference of 
masses of kaons and D mesons. 

IV. REACTION pp -► MM* 



The reactions pp — > .K\ff * and pp — ► DD* are similar to the reactions with KK and DD 
final states, and the corresponding amplitudes are described by the diagrams depicted in 
Fig. 9, where one of the outgoing pseudoscalar mesons M is replaced by the vector one, M*, 
i.e. K — » K*, D — ► .D* etc. Thus, the invariant amplitude for pp — > KK reads 

"ds(*)-| 



rppp^K K* 

1 \j;m i ,n i 



.pp^KK i ,,9k*na9kna w 1 ,,» 



c'(^k;,(m) : 



, (37) 



SO ^ V -Spp-.KK* . 

where A/ is the polarization of the outgoing K* and the other notations have been introduced 
already in Sects. II and III. The baryon exchange trajectories are the same as in the previous 
Sect. Ill, Sp p .j^ K * = Sp p: K K , and C'(t) is defined in Eq. (35) 
The spin dependent amplitude M. has the following form 



M 



PP^rKK 

\t\m 



V 



;>,0=^(M)rLm,r 



(38) 
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with 

^x rm = Vn, [75 (fr - M Y ) (y + 2{ M™+M Y ) {rfi>K - " ^ 7 " )} ] W ^ ' (39) 

where e^. is the polarization vector of the K* meson and 



with 



\f( s t ) = F °°^ 



F 2 (s,t) = Tr[r,r„t] (-<,»» + && 



2sM 2 
F ™ = -Tjf ( M n + M V + 6M N M v z + M 2 (2M 2 N + M 2 )z 2 ) , (40) 

My 

where My = Mk*, Pv = Pk* and z = kvny/(Mn + My). 

In the reaction pp — > K~K* + the total amplitude is the coherent sum of A and E° 
exchange channels. In the case of K°K*°, the amplitude is defined by the S + exchange 
trajectory. 

The amplitude for the K*K reaction has a similar form: 



r M = v 

\f,rrnrn n i 



« + %0fW) ^ - **"™ tf- ' *fr> 75 



« m ,C ■ (41) 



The generalization for DD* and D*D final states may be done in a straightforward 
manner, similarly to the previous Sections. 

1. Differential cross sections 

The differential cross sections of the reactions pp — > KK* and pp — >■ K*K) are exhibited 
in Fig. 15. The left panel shows our prediction for the differential cross sections of the 
reactions pp — > K~K* + , and pp — > K°K*° as a function of t max — £ at initial momentum 
Pl = 10 GeV/c. The dependence of the differential cross sections on energy y/s at fixed 
tmax — t = 0.2 GeV 2 is presented in the right panel. At large energies the cross sections 
behaves similarly to the cross sections of the pp — > KK reactions. 

Our prediction for the differential cross sections of the DD* pair production is presented 
in Fig. 16. The left panel illustrates the dependence on t max — t at fixed pi = 15 GeV/c. The 
right panel exhibits the dependence on the energy excess As 1 / 2 at fixed t max — t = 0.2 GeV 2 . 

20 




0.0 0.2 0.4 0.6 0.8 



t max -t (GeV^) 



10 J 
10 2 

| 10° 

I10" 1 



10" 



10" 



pp->KK t max -t=0.2 GeV : 




2 4 6 8 10 12 
s 1 ' 2 (GeV) 



FIG. 15: Left panel: The differential cross sections of the pp — ► K~K* + (solid curve) and pp — > 
K°K> (dashed curve) reactions as a function of £ max — £ at pl = 10 GeV/c. Right panel: The 
differential cross sections as a function of the energy y^ at £ max — £ = 0.2 GeV 2 . 
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FIG. 16: Left panel: The differential cross sections of the pp — ► D°D* (solid curve) and 
pp — > D~D* + (dashed curve) reactions as a function of £ max -t at pi = 15 GeV/c. Right panel: 
Differential cross section as a function of the energy excess As 1 ' 2 at £ max — t = 0.2 GeV 2 . 

The ratio of the cross sections with D~D* + and D°D final states is defined by the coupling 
constants in KNTi, K*NH and KNA, K*NA interactions and is close to 0.03. The cross 
sections decrease with energy similarly to the pp —+ DD reactions, and therefore, the region 
with small excess energy is more suitable for studying these reactions. 

2. Longitudinal asymmetry 
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Let us consider first the case of the forward production angle (or t = t max ), in pp — » KK* 
(DD*) reactions with pure vector coupling in V*NH (V = K*, D*). In case of large energies, 
where pv — P P , the amplitude has the form 



Txi-mum ~ R(s) (A8 mini + B 6-rmrn) S. 



\i\\. 



(42) 



where Ay is the polarization of the outgoing vector meson, Aj = raj + Hi, and 

M N 



A~V2, 



B 



Mi 



v 



(43) 



where My is the mass of the vector meson. This results in 

M 2 N - 2M V 
M 2 N + 2M y ' 



(44) 



Thus, for the KK* final state, where My ~ Mjv the asymmetry has values A — —0.3 and 
increases at finite production angles because of the spin-orbital interaction. In case of the 
pp — > DD*) reaction, the asymmetry is much smaller: A ~ —0.8 and it also increases with 
the production angle. The finite tensor coupling changes these predictions, especially in the 
charm sector with large My. 
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FIG. 17: The longitudinal asymmetry for the reactions pp —> K~K* + (solid curves) and pp — ► 
K°K*°, (dashed curves). The dot-dashed curves correspond to the calculation with zero tensor 
coupling constant. Left panel: The asymmetry as a function of momentum transfer t max — t at 
PL = 10 GeV. Right panel: The asymmetry as a function of the energy y^ at i max — t = 0.2 GeV 2 . 
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In Fig. 17 we show our prediction for the pp — > KK* reactions. The left panel exhibits 
the t dependence at pl = 10 GeV/c, whereas the right panel depicts the y/s dependence 
at £ max — t = 0.2 GeV 2 . One can see an increase of the asymmetry with —t and its almost 
constant value at large y/s and fixed £ max — £. The difference in A for K°K *° and K + K*~ final 
states is mainly due to the difference in tensor couplings in K*NH and K*NA interactions. 
For completeness, we also show the result for a calculation without tensor couplings. At 
zero production angle the asymmetry is close to our above qualitative estimate. 
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FIG. 18: The longitudinal asymmetry for reactions the pp — ► D°D*° (solid curves) and pp — ► 
D~D + (dashed curves). The dot-dashed curves correspond to the calculation with zero tensor 
coupling constant. Left panel: The asymmetry as a function of momentum transfer £ max — £ at 
PL = 15 GeV. Right panel: The asymmetry as a function of the excess energy A.y/s at £ max — t = 
0.2 GeV 2 

The longitudinal asymmetries for pp — > DD* are presented in Fig. 18. The left panel 
displays the t dependence at initial momentum pl = 15 GeV/c. The right panel shows 
the dependence on the energy excess at t max — t — 0.2 GeV 2 . We also show results for a 
calculation with pure vector couplings. At zero production angle the asymmetry (left panel) 
coincides with our qualitative estimate (cf. Eq. (44)). The effect of the tensor interaction is 
rather large. One can see a big difference between D°D*° and D~D* + final states because 
of the difference in the corresponding tensor couplings. Numerically, result for the D~D* + 
final state is close to a calculation with zero tensor coupling. 
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V. SUMMARY 

In summary, we have analyzed the open charm production in the exclusive binary reac- 
tions pp — * Y C Y C , pp — ► DD and pp — ► /}/)* at small momentum transfer. Our considera- 
tion is based on a modified Regge type model, motivated by quark-gluon string dynamics. 
The most important parameters of the model are the effective charmed meson and baryon 
exchange trajectories and the energy scale parameters. They are found from a consistent 
approach based on the topological decomposition and factorization of the corresponding pla- 
nar quark diagrams. The coupling constants are taken to be the same as in corresponding 
strangeness production reactions, assuming SU(4) symmetry. Unknown residual functions 
are found from the comparison of the calculation of pp — > AA (AS) and pp — > KK reactions 
with available (although old) experimental data. As a result, we obtained the absolute value 
of the corresponding cross sections in the energy range of future FAIR experiments. 

For the first time we made predictions for the longitudinal asymmetry, which is quite dif- 
ferent in different processes with non-trivial t and s dependencies. In each case we presented 
an analytical estimate for the forward production with the aim to understand the physics 
of the asymmetry. 

Our calculations of cross sections and longitudinal asymmetries in the exclusive reactions 
pp — ► YY, pp — ■» KK and pp — * KK* have also an independent interest for forthcoming 
experiments at FAIR as a first prediction of the open strangeness production in peripheral 
reactions in this energy region. Our consideration may serve as a first step towards more 
involved reaction mechanisms and extensions to pp collisions. 
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